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ABSTRACT 

Using constrained simulations of the local Universe for generic cold dark matter and for 1 keV 
warm dark matter, we investigate the difference in the abundance of dark matter halos in the local 
environment. Wc find that the mass function within 20 ft. _1 Mpc of the Local Group is ~ 2 times larger 
than the universal mass function in the 10 9 — 10 13 /i _1 M© mass range. Imposing the field of view 
of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function 
in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, 
employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare 
the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find 
that our simulations are able to reproduce the velocity function in the 80 — 300 km s -1 velocity range, 
having a value ~ 10 times larger than the universal velocity function in the Virgo-direction region. 
In the low velocity regime, 35 -80 km s _1 , the warm dark matter simulation reproduces the observed 
flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a 
steep rise in the velocity function towards lower velocities; for V rna x — 35 kms , it forecasts ~ 10 
times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results 
indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions 
about energetic feedback in dwarf galaxies. 

Subject headings: Cosmology: dark matter, galaxies: Local Group 



1. INTRODUCTION 

Despite the success of the ACDM paradigm in describ- 
ing the large scale structure of the Universe, as has been 
attested by improved measurements of the CMB and 
the large scale cluster ing of galaxies ([Seliak et alJl2006bt 
iKomatsu et all 120081) . some potential problems remain 
for the model at smaller scales. One of these challenges 
is related to the prediction of the abundance of low mass 
galaxies. Numerical simulations within the ACDM model 
predict many more dark matter subhalos inside galactic- 
sized host halos than the actual numb er of observed satel- 
lite galaxies aroun d the Milky Way (iKlypin et aTlll 999; 
Moore et all 119991 : iDiemand et"afll2007t ISpringel et all 
20081 ). This discrepancy has been known as the "miss- 
ing satellite problem" . Although some astrophysical pro- 
cesses can be responsible for the suppression of galaxy 
formation in sm all mass halos and lead to the solutio n 
of the problem (ISomervillel[200l iBenson et al.ll2002bl lal: 
iGnedin fc Kravtsov 2006; Koposov et al.ll2009| ). it is also 
possible to alleviate the discrepancy by adopting alter- 
native models with a different type of dark matter. 

Warm dark matter (WDM) models, for instance, 
predict sign ifican tly less su bstructure within halos 
i|Coh'n et alJl2000HBode et al.ll2001l ). The free streaming 
length for dark matter particles depends on their intrinsic 
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properties and establishes a cut-off scale for halo masses 
below which primordial density perturbations are wiped 
out. In the case of neutralinos, one of the favorite candi- 
dates for CDM m odels, with m x ~ lOOGeV, the scale is 
roughly 10" 6 M (jDiemand et alll2005fh For W DM can- 
didat es, such as gravitinos with mg ~ IKeV (see[Stcffcn 
120061 for a recent review on gravitino cosmology), this 
scale is roughly 10 10 M Q (see section 2). Therefore, using 
the same cosmological parameters but different cut-off 
scales results in different predictions for the halo mass 
function (MF). At the low mass end the AWDM mass 
function is expected to be significantly lower than the 
ACDM mass func tion. Above the cut-off scale the differ - 
ences vanish (e.g. lBode etaLll2001tlBarkana et al.ll2001f) . 

Since for high density environments both models pro- 
vide, conceptually different, solutions for the missing 
satellite problem it is also important to test their pre- 
dictions for isolated systems in low density environments 
which are not affected by processes such as ram-pressure 
or tidal striping. A comparison between AWDM and 
ACDM is then free from these astrophysical phenomena 
and conclusions drawn from it are more closely related 
to the nature of dark matter itse lf. Such an approac h 
has been employed for example bv lBlanton et alJ {2008). 

Even though low mass galaxies are difficult to de- 
tect, surveys such as the o ngoing Arecibo Legacy Fast 
ALFA (ALFALFA) survey (|Giovanelli et alJ l2005al ) are 
promising to make such comparison. ALFALFA is ex- 
ploring the local HI Universe with a sky coverage of 
7000 deg 2 and it is designed to detect objects with 
masses as low as 3 x 10 7 Mf7> at a Virgo- cluster dis- 
tance (11.7 h^Mpc) (|Giovanelli et al.ll2007D . The sur- 
vey is divided into two regions on the sky. One within 
the solid angle 07 h 30 m <R.A.(J2000.0)< 16 ,l 30 m and 
0° <DEC(J2000.0)< +36°, which includes the Virgo 
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cluster, and another region within the same declination 
range and 22 h <R.A.(J2000.0)< 03\ In the reminder of 
this work, we will loosely refer to these two regions, ad- 
ditionally constrained to distances less than 20 /i _1 Mpc, 
as the "Virgo-direction region (VdR)" and "anti- Virgo- 
direction region (aVdR)" respectively. We note that the 
catalogs released by ALFALFA so far in these two re- 
gions effectively correspond to a high and a low density 
environment. The HI sources detected by ALFALFA are 
typically associated with star forming disk galaxies, spi- 
rals in rich clusters are less likely to be detected sinc e 
these galaxies might be HI deficient (jSolanes et alf 2001). 

In order to make predictions on the abundance of 
low mass galaxies based on the ACDM and AWDM 
models for our local environment, it is advantageous 
to use co nstrained simulations (CSs) of the local uni- 
verse fe.g.lBistolas fc Hoffman! [l99l iMathis et al.ll2002t 
iKlypin et al.l I2003TK These simulations are constructed 
to reproduce the gross features of the nearby Universe, 
such as the Local Supercluster and the Virgo cluster. 
This can be achieved by setting the initial conditions of 
the simulations as constrained realizations of Gaussian 
fields, where actual observational data are used to impose 
the constraints. By using such simulations, biases which 
are present because of the particular structure of the lo- 
cal environment, are minimized and a better comparison 
between simulations and observations can be achieved. 

Although the abundance of dark matter halos is usu- 
ally quantified by the mass function, which describes the 
number density of halos as a function of mass, usually 
presented either in a differential o r integral form (e.g. 
iReed et all l2007t iLukic et aLll2007t ). an alternative ap- 
proach is to use the velocity function, similarly defined 
as the number density of halos as a funct i on of maxi- 
mum circular velocity (e.g. iGonzalez et al.l [20001. The 
advantage of the velocity function (VF) is that it can be 
compared more directly with observational results since 
it avoids the more complicated problem of relating dark 
matter halo masses to galaxy luminosities. Instead, the 
velocity function of galaxies can be theoretically esti- 
mated using the velocity function of halos by incorpo- 
rating a model of baryonic infall; the only processes that 
affect the velocity function are those that modify the 
gravitational potential of galactic systems. In the present 
work we show analysis of both, the mass and velocity 
functions, but concentrate on the latter for comparisons 
with the ALFALFA survey. 

The objective of our paper is to run a ACDM and a 
AWDM simulation of the local environment where the 
same set of constraints have been imposed and to com- 
pare the abundance and distribution of dark matter ha- 
los for the two cosmologies. In addition, the results are 
used to predict the abundance of HI sources that will be 
detected by the on-going HI blind survey ALFALFA. 

The paper is organized as follows. In section 2 we de- 
scribe the setting of the simulations. The definition of 
the coordinate system that we use to impose the field 
of view of ALFALFA in the CSs is described in section 
3. In section 4 we present results on the abundance of 
dark matter halos. In section 5 we analyze in particu- 
lar the abundance for the restricted Virgo-direction and 
anti- Virgo-direction regions cataloged by the ALFALFA 
survey to date and present predictions on the velocity 
function of HI sources. The summary and conclusions of 



our work are given in section 6. 

2. CONSTRAINED SIMULATIONS OF THE LOCAL 
ENVIRONMENT 

We chose the cosmological parameters for our sim- 
ulatio ns to be cons i stent with the WMAP 3-year re- 
sults (|Spergel et alJ 120071 ): fi m = 0.24, n A = 0.76, 
H = lOO/ikms^Mpc" 1 with h = 0.73, n = 0.95 and 
as = 0.75. For both cosmologies the theoretical (uncon- 
strained) linear power spectrum at z = is shown in Fig. 
Q] (blue for CDM and red for WDM). The CDM power 
spectrum was computed from a Boltzmann code by W. 
Hu and was kindly provided to us. 

2.1. AWDM simulation settings 

The WDM power spectrum was computed by rescal- 
ing the CDM power spectrum using a fitting function 
that approximates the transfer function for a thermal 
WDM particle with mwDM = 1 keV. The filtering scale 
(or free-streaming length) for this WDM particle mass is 
SSOh^lspc, or about - 1.1 x 10 10 /i _1 M e f or the filter- 
ing m ass (following the definitions given by iBode et"al 
120011) . We used the fitting function from IViel et al. 
(|2005h (see t heir eqs. 5- 7 ) whi ch is very similar to the 
one given by IBode et al.l J2001), however, according to 
the authors, it is a better fitting formula for Boltzmann 
codes. 

As stated by recent Lyman-a forest, CMB and ga laxy 
clus tering data (IViel et al.ll2 006; Scli ak et"aLll2006al ). see 
also iMiranda fc Macciol (|2007l K the lower limit for the 
mass of thermal WDM candidates, for the case where 
WDM is the dominant form of dark matter, has been 
reported to be mwDM ~ 2keV at the 2a level. However, 
these estimates for thermal relics may b e contaminated 
by syst ematic errors (see for example iBovarskv et al.l 
(|2008al) for a discussion on the complications associ- 
ated to the Lyman- ct forest me t hod). For instance, the 
constraint given in IViel et all ((2006) would reduce to 
ttiwdm > 0.9 keV if the highest redshift bins of the 
Lyman-a forest data are rejected from the analysis and 
only the more reliable data ba sed on z < 3.2 i s taken 
into account. In a recent paper, IBovarskv et~aTl (|2008bD 
revisit the lower bounds on the mass of WDM particles 
and find my/DM > 1.7keV (2a level). However, after dis- 
cussing with detail the systematic uncertainties in their 
method they conclude that the mass bounds are reliable 
within ~ 30% of uncertainty. All these analysis put our 
choice of raw dm = 1 keV close to the most recent lower 
bound, but it is a choice that is still not ruled out. 

The velocity dispersion of the WDM particles needs, in 
principle, to be introduced in the initial conditions for the 
simulations since after all, it is what causes the smearing 
of small scale primordial perturbations. This is usually 
done by introducing a random veloci ty field according 
to a Fermi-Dirac distribution function .IBode et al.ll200lT 
IColm et al.ll2008T ). However, if this velocity dispersion 
field is introduced at random, a certain degree of white 
noise is generated (shot noise) due to the finite number of 
simulation particles which produces spuriou s small scale 
powe r in the WDM spectrum (see Fig.l of IColm et al.l 
2008). The amplitude of the rms velocity of the random 
component to be added depends on the nature of the 
WDM particle and on the redshift of the initial condi- 
tions. For thermal relics, it is larger for higher redshifts 
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and smaller particle masses. For z = 50 the rms velocity 
is ~ 2.2 kms " 1 for mwDM = 1 keV (following the for- 
mula given bv lBode et al.|[200ll ). far lower than the veloc- 
ities induced by gravitational collapse of structures hav- 
ing scales larger than the minimum scale we can resolve 
within our simulations. Therefore we do not incorporate 
any random velocity field into our initial conditions. 

To back our approach we computed the comoving 
Jeans length associated with a velocity dispersion of 
2.2 kms -1 . We found that the time to build up pres- 
sure support against gravitational collapse is similar to 
the collapse time for comoving scales of ~ 33 h~ 1 kpc, 
corresponding to a Jeans mass of ~ 10 6 fc _1 Mg. Differ- 
ences in the structure formation are only expected for 
masses of the order and below the Jeans mass, which is 
lower than the particle mass in our simulations (see be- 
low). Thus, the effect of thermal velocities is too small 
to have a notable effect on the abundance of low mass 
halos in our simulations. St ill, they could hav e an impact 
in the inner part of halos (|Colm et alj 120081 ). However, 
since the shot noise produced by their introduction has a 
large spurious effect, we decided not to include them in 
our simulations. For the purposes of this work, this has 
no consequences in our analysis. 

2.2. Constrained simulations 

The initial conditions for the constrained simu lations 
were set up using the (Hoffm an fc Riba k 1991) algo- 
rithm of constrained realizations of Gaussian random 
fields. Two types of data sets are used as input for 
the algorithm. The first data set is made of radial ve- 
loci ties of galaxies dra wn from the catalogs: MARK 
III (IWillick et al1IT9 97). surface brightness fluctuations 
dTonrv et al.l 2001) and the C atalog of Nearby Galax- 
ies ( Karachentsev et al.l |2004) . Present epoch peculiar 
velocities are less affected by non-linear effects and are 
therefore imposed as linear constraints on the primordial 
perturbation field (|Zaroubi et al.lll999D. Thi s appr oach 
follows the CSs per formed by lKravtsov et ail (|2002D and 
iKlypin et alj ((2003) . The other data set is obtained from 
the catalog of nearby X-ray se lected clusters of galax- 
ies ([Reiprich fc Bohringerl 12002). Assuming the spheri- 
cal top-hat model and using the virial parameters of a 
cluster the linear over-density of the cluster is derived. 
The estimated linear overdensity is imposed on the virial 
mass scale of the cluster as a constraint. The density 
and velocity fields on scales larger than 5/i _1 Mpc are 
strongly constrained by the imposed data. 

Using the above initial conditions, we carried out th e 
simulations using the code GADGET-2 (|Springelll20"05T ). 
Both simulations follow the evolution of 1024 3 dark mat- 
ter particles from z = 50 to z — in a box of size 
L = 64 /i -1 Mpc. The associated Nyquist frequency and 
fundamental mode are represented with vertical lines in 
Fig. [U The dark matter particle mass is tudm — 
1.63 x 10 7 /i _1 Mq. The simulations were started with 
a fixed comoving softening length (Plummer- equivalent) 
of e = 1.6 /i _1 kpc. Once the corresponding physical co- 
moving softening reached e = 0.8 /i _1 kpc, it was kept 
constant at this value. 

Fig. [2] displays the projected dark matter distribution 
of the CSs at z = for the ACDM and the A WDM 
models respectively. The highest and lowest concentra- 
tions are colored red and black respectively. Both pan- 



els show the matter distribution within a slice that is 
8 /i -1 Mpc thick projected onto the X-Y plane and cen- 
tered on Z — 24 ft, -1 Mpc. The figure captures the Local 
Supercluster (LSC) which is the filamentary structure 
crossing the image plane horizontally. It is the most 
prominent feature. The locations of the virtual Local 
Group (LG) 6 and the Virgo cluster are marked as well. 
A description of the identification of these objects will fol- 
low below. A visual inspection of the two figures already 
reveals a deficit of small scale structure in the AWDM 
simulation, as expected. 

Halos in the simula tions were identified using AHF 7 , 
AMIGA'S halo finder (jKnollmann fc Knebdl2009f) . AHF 
identifies halos as local density maxima in an adaptively 
smoothed density field using a hierarchy of grids and a 
refinement criterion. The latter was chosen so that a 
grid cell is subdivided until each subsection contains less 
than 5 particles. In this way, the size of the smallest cells 
is comparable to the force resolution of our simulations. 
Halos, and/or subhalos, are formed by groups of par- 
ticles that are gravitationally bound to a given density 
peak. For each halo in the final catalog, AHF provides 
a list of internal properties, the most important quan- 
tities for the current study are: the viral radius r« r , 
defined as the radius that contains a mean overdensity 
p{fvir) = Ap cr it, where p cr u is the critical density and 
we have chosen A ~ 94 at z — 0, a value computed 
according to the spherical collapse model for our cosmo- 
logical parameters 8 ; the corresponding virial mass M V i r ; 
the maximum rotational velocity V max ] and the spin pa- 
rameter A = J^JW\ I GM„£ (iPeeblesI 119691) . where J 
and E are the total angular momentum and energy of the 
halo. Subhalos identified within each halo were removed 
from the catalog; for the reminder of our analysis we will 
use main halos only. 

3. DEFINITION OF THE COORDINATE SYSTEM 

The significance of our study relies on a proper sim- 
ulation of the local environment. However, the con- 
straints imposed on the initial conditions of the simu- 
lations leave some freedom for the evolution, mainly on 
small scales (< 5 /i -1 Mpc). This can partly be accounted 
for by a careful adjustment of the coordinate system. We 
have done so using the following steps. First we identi- 
fied an appropriate LG within the CSs. To that pur- 
pose w e followed the crit eria described in iMacciq et all 
(|2005f) (see also table 2 of lMartinez-Vaauero et al.ll2007t) 
to identify LG candidates within the CSs. These are i) 
match in mass, proximity and kinematics of a MW-M31- 
halo-like binary system, i.e., we look for a pair of ha- 
los with a maximum rotational velocity, V max , between 
125 kms -1 and 270 kms -1 , with a distance between the 
pair < 1 /i _1 Mpc and with a negative relative velocity; 
ii) absence of a massive near-by halo, i.e., no halos with 
masses larger than the members of the pair within a ra- 
dius of 2 /i _1 Mpc; and iii) presence of a Virgo-like halo, 
500 < Vmax < 1500 kms -1 , at the appropriate distance, 
5 - 12 /i- x Mpc. 

In addition, we favor LG candidates which are located 

6 In what follows, the LG refers to a pair of halos associated 
with the MW and M31 galaxies 

7 http://www.aip.de/People/AKnebe/AMIGA/ 

8 For Q m + Q A = 1, A sa 178Q^ 45 fe.e. lEke et al||l998h 



4 



Zavala et al. 



close to the center of the simulation box since the con- 
strained initial conditions place the LG progenitor ex- 
actly at the center. Subsequent dynamical evolution may 
displace the entire environment by some Mpcs. 

Slices containing the best LG candidates for the CDM 
and WDM run are shown in Fig. [2] where we have marked 
the location of the LG and the Virgo cluster. If not stated 
otherwise, "LG" and "Virgo" denote these best possi- 
ble candidates. Table 1 contains some properties of the 
main objects in our simulations: the LG, and the clus- 
ters Virgo and Fornax (see below). Observed estimates 
for these properties are also given in the table. 

Based on the locations of the LG, Virgo and the local 
environment we now aim to define a superg alactic co- 
ordinate system rtde Vaucouleurs et al.l (|1991f k see also 
iLahav et al.ll2000h . We will use such coordinate system 
as the basis to identify the ALFALFA regions in the (sim- 
ulated) sky. To define it, we first assume that the equa- 
torial plane of the supergalactic coordinate system lies 
in the supergalactic plane (SGP), which is spanned by 
the LG and the LSC. Thus, besides Virgo, we need to 
find another cluster belonging to the LSC to mathemat- 
ically define the LSC. As revealed by Fig. [2j there is a 
prominent cluster on the right hand side of Virgo. Its 
location and mass closely resemble those of the observed 
Ursa Major cluster (1.6 x 10 13 /i -1 M Q at a distance of 
~ 11.1 h Mpc from the LG). Albeit this choice seems 
natural we have also tested other clusters within the sim- 
ulated LSC to define alternative coordinate systems. The 
final results turned out to be almost independent of this 
choice. Using this SGP we construct an orthonormal vec- 
tor basis with the origin placed at the MW and rotate it 
until the simulated Virgo cluster is located at the same 
longitude of the observed one, i.e., at a supergalactic 
longitude (SGL) of 102.45°. Finally, we tilt the plane to 
achieve a supergalactic latitude (SGB) of 2.84° to match 
the one that is observed. With this procedure, the sim- 
ulated Virgo is located at the same angular position as 
the observed one. Below we will introduce an alternative 
coordinate system. To avoid confusion we refer to this 
first one as SG zero . 

The result of our coordinate definition is visualized in 
the left panel of Fig.[3j It shows a sky map with the angu- 
lar distribution of halos in Equatorial coordinates (RA, 
DEC) in a Mollweide projection for the ACDM simu- 
lation. Only halos within 20/i _1 Mpc from the origin 
(simulated MW) are included. The size of the sphere is 
limited by the size of the simulation box. Larger radii are 
likely to contain spurious structures caused by periodic 
boundary conditions. The map was created using the 
HEALPIX software 9 using N pix = 12 (64) 2 pixels. The 
value of each pixel is given by a mass-weighted count of 
all halos located in that pixel. Afterwards, a smooth- 
ing of the map was done using a Gaussian beam with a 
FHWM of 7°; the different color scale is a measure of 
the values in the map, from red-to-blue for high-to-low 
values. The angular positions of all halos with masses 
larger than 5 x 10 9 /i -1 M Q also appear in the map as 
black points. The voids and high density regions within 
the local simulated volume can be clearly appreciated in 
the sky map. The prominent high density region cross- 
ing the whole map almost vertically at the center is the 

9 http://hcalpix.jpl.nasa.gov/ 



simulated LSC with Virgo roughly in the middle (black 
circle). By construction, the location of the real Virgo 
(black square) is identical with the simulated one. The 
boxes in the center and on the sides of the map give the 
boundaries of the VdR and aVdR, respectively, accessi- 
ble to the ALFALFA survey. 

Since the LSC is roughly in place in our CSs and be- 
cause the simulated Virgo and the real one are at the 
same angular position in the sky (and almost at the same 
distance from the LG) , we believe that the procedure de- 
scribed in the last paragraph led to an appropriate co- 
ordinate system allowing us to use our CSs to simulate 
the VdR of the sky surveyed by ALFALFA. However, a 
visual impression of the aVdR (enclosed by the boxes on 
the sides of the sky map in the left panel of Fig. [3]) indi- 
cates tentative problems with this simulated area of the 
sky. It contains a significant part of a filamentary struc- 
ture with a high density of halos, stretching from around 
(0°,0°) to (45°, +20°). Such structure is associated with 
the simulated Fornax cluster. This cluster was imposed 
as a part of the simulation constraints, see Table 1, it 
appears as a black circle in the right corner of the sky 
map. Its angular position is however out of place, the 
real angular position of the Fornax cluster is marked as 
the lower right black square in Fig. [3] 

Such deviation is within the expected variations of CSs 
due to intrinsic uncertainties. For instance, the velocity 
constraints used to produce the CSs have still large er- 
rors, making the random component of the CSs more sig- 
nificant. Also, the present day positions of the massive 
clusters obtained from observational data are imposed 
on the initial conditions of the CSs. Both of these con- 
ditions imply that the dynamical evolution of the sim- 
ulations shifts the positions of the clusters and modify 
the structures finally obtained in the CSs. Furthermore, 
since scales smaller than ~ 5 ft. Mpc are unconstrained, 
the position of the LG is not imposed directly by the 
constraints. Finally, to be able to explore low mass halos 
we have used a small box size for the CSs; since the LSC 
stretches from one side of the box to the other, periodic 
boundary conditions distort the shape of the LSC. De- 
spite all these difficulties, we arc confident that the CSs 
and the choice we have made for the location of the LG 
are reliable enough to make the comparison we intend in 
this work. 

The result to keep in mind is that Fornax appears in- 
side the aVdR whereas in reality it does not. Then, by 
using SG zero we would over-predict the abundance of ha- 
los in this region since by mistake our simulated survey 
probes a region of higher density than the one expected 
from observations. 

To ameliorate this problem, we carried out an addi- 
tional adjustment of the coordinate system. We relax 
the requirement that the angular position of the sim- 
ulated Virgo has to coincide exactly with the angular 
position of the real Virgo. Instead, we look for a coor- 
dinate system with the same origin, and that minimizes 
the quadratic sum of the distances between the simulated 
and real clusters Virgo and Fornax. We refer to this coor- 
dinate system as SG m i n . The distribution of halos after 
this rotation can be seen in the sky map at the right 
panel of Fig. [3] In the coordinate system SG m i n , the 
distances between the real and simulated Virgo and For- 
nax clusters are 4.5 ft _1 Mpc and 5 /i _1 Mpc respectively. 
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The main effect of this rotation is that the (ALFALFA) 
aVdR region no longer includes significant parts of the 
filamentary structure associated to the Fornax cluster. 

We are confident that the adoption of the new coor- 
dinate system, SG m i n , is consistent with the freedom 
inherent to the CSs. In section 15.11 we investigate the 
impact of the rotation in a more general sense. There we 
can show that minor rotations, as the one adopted for 
the final adjustment of the coordinate system, change 
the halo abundance in the VdR by less than 30%. Thus, 
the adjustment of the coordinate system introduces only 
a small uncertainty that we will however keep in mind 
for the interpretation of our results. Here we have only 
discussed the adjustment of the ACDM coordinate sys- 
tem, we repeated the same procedure for the AWDM 
simulation. 

4. HALO ABUNDANCE 
4.1. Global mass function 

The blue and red solid lines in the upper panel of 
Fig. [4] display the differential mass function for the 
ACDM and AWDM CSs, respectively. Analytical pre- 
dictions for both cosmologies are shown as dashed 
lines followin g the Sheth fe Torm e n formalism (S-T 
formalism, ISheth fc Toriner] [l99l ISheth et all 12001 : 
iSheth fc Tormenl 12002). For th eir computa t ion w e used 
the public code described in iReed et ail ((2007). In 
this prediction, the only difference between both cos- 
mologies is the suppression of the power spectrum at 
small scales in the AWDM case. The statistical error 
bars presented in the figure are Poisson errors , employ- 
ing the definition given in iLukic et al.l (|2007h : a± = 
^/N + 1/4 ±1/2, where N is the number of halos per 
bin. The value of the filtering mass for the WDM simu- 
lation, 1.1 x 10 10 A^'Mq, is marked in the figure with a 
vertical solid line. 

At the high mass end there is an excess of massive ha- 
los in our simulations. This is caused by the constraints 
which enforce the growth of a very massive structure, the 
LSC, within the relatively small volume of the simulation 
box. At the low mass end we see di screteness features as 
described in IWang fc Whitel (|2007l ). According to their 
study, the limiting mass that can be trusted is given by 
the formula: Mu m = 10.lpdk~^ ak , where d is the inter- 
particle separation, k pea k is the wave number for which 
k 3 P(k) reaches its maximum and p = QdmPctU- For the 
case of our AWDM simulation: M Um = 3 x 10 9 h^M®. 
As can be seen in Fig. 21 this value (indicated by the 
dotted vertical line) marks the mass limit below which 
there is an artificial rise in the WDM MF, an indication 
for the onset of discreteness effects. 

The difference between the MF of the different sim- 
ulations becomes even more clear in the lower panel of 
Fig. [U where we plot the ratio of the measured MFs 
to the value of the S-T MF for the ACDM cosmology. 
Clearly, the abundance of low mass halos in the AWDM 
case (for masses larger than Mn m ) is considerably lower 
than the predictions using the S-T fo rmalism. This dis - 
crepancy has been found before (e.g. jBode et aLlfeOOlf ). 
so one should not expect good agreement between the 
S-T approach in the range of masses close to and below 
the AWDM filtering mass. In the following we will use 
the ACDM S-T mass function as a reference to present 



some of our results. 

4.2. Global velocity junction 

As was mentioned in the introduction, a more di- 
rect comparison for the abundance of structure between 
our CSs and observations of the local Universe can be 
achieved by constructing the velocity function of dark 
matter halos. The maximum rotational velocities are 
used instead of the masses to calculate the abundance of 
halos per logarithmic V max bin. Provided that a physical 
connection between V max and the measured maximum 
rotational velocity of spiral galaxies can be established, 
our velocity function can be directly compared with ob- 
servations of galactic discs. Further below, in 15.21 we 
describe a simplified model which is designed to accom- 
plish this goal. 

The upper panel of Fig. [5] shows the differential VF for 
halos in both simulations (the line styles and colors are 
the same as in Fig. 2]). Also shown is a prediction for 
the VF in the ACDM cas e (dashed lin e ) obta ined using 
the procedure outlined in lSigad et all (|200Ct ). In brief, 
the procedure is the following: i) a "virtual" sample of 
halos is generated in such a way that its MF mimics 
the one given by the S-T formalism; ii) a concentrati on 
value taken from a log-normal distribution (Jing 2 0001 ) is 
assigned to each halo; mean and standard deviation val- 
ues for the log- normal d i stribu tion were taken from the 
analysis bv lMaccio et al.l (|2008l ) for the set of simulations 
with WMAP3 cosmology: (logc) = 1.775 - 0.0881ogM 200 
and aiogc = 0.132; iii) assuming a NFW profile, we com- 
pute V m ax for each halo using the concentration and 
viri al velocity c o rrespo nding to that halo (e.g., see eq. 
7 of ISigad et al . (2000)). A prediction in the case of the 
AWDM simulation is not given since the previously de- 
scribed steps are all uncertain in this case. 

The vertical solid and dotted lines in Fig. [5] are the es- 
timated values for the maximum velocities corresponding 
to the filtering and limiting masses for the AWDM sim- 
ulation. For the computation of the corresponding virial 
velocities we used the relation between the virial mass 
of the halo and its virial radius: M„ r = ^TrAp cr i t r^ ir . 
Then, the virial velocity is simply given by, V^ ir = 
GM V i r /r V i r . Using the same mass-concentration rela- 
tion described in the paragraph above, we compute the 
maximum circular velocities for the filtering and limiting 
mass and find the values of 36 kms -1 and 24 kms -1 . 

The lower panel of Fig. [5] shows the ratio of the mea- 
sured VFs to the analytical one displayed as dashed line 
in the upper panel 10 . This figure is analogous to the ra- 
tios of the MFs as shown in the lower panel of FigJH The 
analytical estimate for the VF follows closely the result 
of the ACDM CS for most of the velocity range. The 
difference at the high mass end is due to the constraints 
imposed in the simulation. 

The downward bend for the ACDM velocity function 
at the low-velocity end is caused by the mass cut-off for 
halos (~ 4x 10 8 /i -1 M Q ). The reason why we do not 
see an abrupt cut off is due to the spread in halo con- 
centrations which causes a similar spread in V max - The 
behavior of the AWDM VF is analogous to the behavior 
of the AWDM MF. Discreteness effects can be seen for 

10 Recall that the analytical result is given by the S-T formalism 
for ACDM only 
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velocities lower than the limiting V max (dotted line) . For 
velocities just above this limiting value, the difference 
between both simulations is approximately an order of 
magnitude. 

4.3. Abundance of halos in the local environment 

We now turn our attention from the global to the lo- 
cal abundance of dark matter halos. In Fig. [5] we show 
the ratio of the differential MFs to the ACDM S-T pre- 
diction for spheres of different size centered on the LG. 
The upper and lower panels show the results for the 
ACDM and AWDM simulations. In both panels the ra- 
tios decrease with increasing radius of the spheres. The 
smallest sphere has a radius of 15 /i _1 Mpc resulting in 
the highest abundance (uppermost curve in each panel). 
The following sequence of curves corresponds to radii in- 
creased by intervals of 5 ft Mpc up to a maximum ra- 
dius of 40 /i _1 Mpc. The curve corresponding to a sphere 
with radius 20 h~ 1 Mpc is highlighted with a thicker dash- 
dotted line. For reference, the result for the MF in the 
whole cubic box is presented as a thick solid line (see 
Fig. Q|. For spheres growing beyond the border of the 
simulation box we use periodic boundary conditions to 
fill in the volume of the sphere. 

According to Fig. [HI the local environment shows an 
overabundance of halos compared to the mean abun- 
dance in the whole box. Within 20 /i _1 Mpc (which was 
the radius used to produce Fig. [3|) the halo abundance 
is about 2 times larger than that in the entire box. In 
what follows all results are restricted to a sphere with 
radius of 20 /i _1 Mpc, which is close to the maximum ra- 
dius that a sphere centered at the LG can have and still 
lie completely inside the simulation box. 

5. PREDICTIONS FOR THE ALFALFA SURVEY 
5.1. Simulated field of view of ALFALFA 

Fig. [7] displays the distributions of halos within the 
ALFALFA field of view projected onto the plane of the 
supergalactic coordinate system, SG m i n , which was in- 
troduced at the end of section [3] The panels on the 
left are for the ACDM simulation and on the right for 
the AWDM simulation. The red and blue dots give the 
position of halos within the VdR and the aVdR respec- 
tively. Only halos with a distance to the LG less than 
20 /i _1 Mpe and with masses larger than 5 x 10 9 h~ 1 M@ 
are shown. The difference in halo abundance between 
both regions is clearly visible. 

Fig. [5] shows the radial dependence of the number den- 
sity of halos, rih, normalized to the total number density 
of halos, n sp hi m the 20/i _1 Mpc sphere. The red and 
blue lines show the result for the VdR and aVdR, and 
the black line for the whole sphere. The upper and lower 
panels arc for the ACDM and AWDM simulations, re- 
spectively. The peak at a distance of ~ 11 /i _1 Mpc in 
the VdR (red line) is caused by halos associated with the 
LSC in the vicinity of the Virgo cluster. For radii larger 
than 10 /i _1 Mpc the VdR is significantly overdense, com- 
pared to the density in the sphere, whereas the aVdR is 
underdense at all radii. 

The difference between VdR and aVdR is also clear in 
Fig. [9] where we show the differential VF for both simula- 
tions (solid lines for ACDM and dotted lines for AWDM) 
in different regions around the LG: the whole 20 h~ 1 Mpc 



sphere (black solid and dotted lines), the VdR (red lines) 
and the aVdR (blue lines). The values for the filtering 
(vertical solid line) and limiting (vertical dotted line) ve- 
locities for the AWDM simulation are also shown for ref- 
erence. 

In both simulations, the difference between the VdR 
and aVdR is clear and qualitatively as expected, the 
former being an overdense region and the latter an un- 
derdense region, relative to the whole region contained 
within the sphere of 20/i _1 Mpc radius. We note how- 
ever, that the aVdR has actually a similar, although 
smaller, density than the mean cosmological density, es- 
pecially for the ACDM simulation. This becomes appar- 
ent from the comparison between the blue solid and black 
dashed lines representing the aVdR and the S-T predic- 
tion, respectively, in the ACDM case (see also Fig. [5J. 
This fact is partially related to the limitations of our CSs 
(already mentioned in detail in section 3). For halo ve- 
locities lower than 100 km s" 1 and larger than the WDM 
limiting velocity (24kms _1 ), the difference between the 
VdR and aVdR is approximately constant in both cos- 
mologies, the ratio of their differential VFs in this range 
is roughly 3. 

To conclude this section we investigate the robustness 
of our results in relation to the choice of the coordinate 
system. For that purpose, Fig. 1101 displays a comparison 
of the VFs for moderate rotations of the supergalactic co- 
ordinate system defined in section [3J The comparison is 
done using the ratio of the velocity functions to the veloc- 
ity function in the whole sphere with 20 /i _1 Mpc radius 
(VF/VF20). The black and blue dashed lines show the 
ratios of the VFs for the VdR and aVdR based on the ini- 
tial coordinate system, SG zero . This coordinate system 
was set up to optimize the agreement between the real 
and simulated Virgo location, but failed to get the For- 
nax cluster in place. In this system, the simulated aVdR 
contains a significant part of the filamentary structure as- 
sociated with the Fornax cluster (see left panel of Fig. [3]) , 
which is not being surveyed in the observed field of view. 
Therefore, the abundance of halos in the simulated aVdR 
is unexpectedly high. The solid black and blue lines show 
the same quantities using SG m i n . With it, the simulated 
aVdR does not contain Fornax any longer and the halo 
abundance decreases substantially. The results based on 
the systems SG zero and SG m in differ by ~ 15 — 20% for 
the VdR and by ~ 40-50% for the aVdR. The red region 
shown in Fig. [TOl displays the ratio of the VFs functions 
in the VdR for rotations of system SG zero up to 20° in 
both supergalactic latitude and longitude. This figure 
indicates that the volume density of halos in the VdR 
stays within 20 — 30% of its value in the coordinate sys- 
tem SG zero for moderate rotations around it and halo 
velocities below 100 kms" 1 (which is the range we are 
ultimately interested in). The velocity function in the 
VdR is therefore robust against moderate rotations of 
the coordinate system due to the appearance of adjacent 
high and low density regions. For the aVdR, the results 
are more sensitive to rotations of the original coordinate 
system. These general results for the choice of coordinate 
system apply in the AWDM case. 

5.2. Velocity function of disk galaxies 

In order to directly compare our simulations with ob- 
servations we need to populate the dark matter halos 
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with galaxies. Essentially, we need to use a method to 
connect the maximum circular velocity measured for disk 
galaxies with the properties of the hosting halo. Since a 
full semi-analytical treatment goes beyond the scope of 
the current study, we use a simplified scheme that, nev- 
ertheless, allows us to make predictions on the velocity 
function of galaxies in the local environment. First, ap- 
propriate halos need to be selected. The high-velocity 
(large-mass) halos are associated with groups and clus- 
ters of galaxies. In our scheme we exclude halos with 
masses larger than 10 13 h~ 1 M Q since we are interested in 
low mass isolated galaxies. Next, we assume each of the 
remaining halos to contain only one disk galaxy (for iso- 
lated galaxies in the local Universe, the fraction of spiral 
galaxies lies between 80-90% , e.g. ISulentic et al.l [20061 : 
iHernandez- Toledo et al.|[2008f ). Although many of these 
halos contain significant substructures, which in princi- 
ple can be populated with satellite galaxies, we restrict 
the present study to main halos and the central galaxies 
within them. This can be justified because the fraction 
of sa tellite galaxies lies between 10 — 40% (|Zheng et al.l 
|2007|) and approximately 50% of them are unlikely to be 
detected in surveys like ALFALFA. The latter estimate 
is roughly correct because, for optically selected sam- 
ples , a large fraction of satellite galaxies have red col- 
ors (IWeinmann et al.ll2006t IWang et alJl20"07l iFont et all 
|2008| ) and have probably already l ost most of their gas . 
HI studies of nearby groups (e.g. iKilborn et al.1 120051 ) 
similarly show that about half of the satellites would lay 
below the sensitivity threshold for HI detection. There- 
fore, our results may underestimate the abundance of HI 
sources by 5 — 20%, and thus, the lack of satellites in our 
modelling is not a major source of uncertainty. 

We compute the circular velocities of the disks lying 
at the center of each halo using the analytical model of 
I Mo et ah! (1998). The principal hypotheses of this model 
are the conservation of specific angular momentum of 
both, the dark matter and the gaseous components, and 
the equality of disk and halo specific angular momenta 
(Jdisk/Mdisk = Jh/M h ) during the process of disk forma- 
tion. These assumptions are key conditions for the for- 
mation of realistic disks in hydrodynamical simulations 
(|Zavala et al.| [2008). The rotation curve of the galactic 
system (disk+halo) is given by the combined gravita- 
tional effects of the disk (which ends up with an expo- 
nential surface density profile) and the adiabatically con- 
tracted halo. The maximum rotational velocity (V max : 
disk+halo) of the disk and the maximum rotational ve- 
locity of the pure dark matter halo (V ma x,h) are related 
by V m ax = G( A, f dls k ) V ma x,h , the function G(X, fdisk) 
depends on the spin parameter of the halo, A, and the 
fraction of baryonic mass that is used to assemble the 
disk, f disk = M d /M vir . The V max /V ma x,h ratio increases 
with fdisk and decreases with A. We note that this ra- 
tio is nearly independent of the halo concentration; for 
fixed values of fdisk and A, variations of c within the typ- 
ical scatter in the mass-concentration relation produce a 
change of less tha n 2%. The fun ction G(A, fdisk) has been 
approximated in IZavalal (|2003f ) with a fitting function 
that has an accuracy larger than 96% for A G [0.02, 0.1]: 



G(X, fdisk) = 1M 1- 



0-ll/d»fc + 5 x 10- 
A 



(1) 



The parameter fdisk actually depends on the processes 



occurring during galaxy formation, such as gas cooling 
and feedback. For simplicity we take a constant value 
fdisk = 0.03 for all galaxies. However, at the end of the 
following section we address the effects of SN feedback. 
For the constant value of fdisk, we take Xu m — 0.02 as a 
limiting value, b elow which the disk would be unstable 
(|Mo et al.lll998[ ). We note however that the large ma- 
jority of the halos in our CSs have A > \u m , halos with 
lower values are not considered in the analysis. 

Applying the scheme described above we obtain the 
velocity function of the modeled disk galaxies. The re- 
sult is shown in Fig. QTJ colors and line styles are the 
same as those used in Fig. [9] We have included in the 
figure the values of the filtering and limiting velocities 
for halos related to the AWDM (vertical solid and dot- 
ted lines). These values mark lower limits for the corre- 
sponding quantities in the case of modeled disks in the 
AWDM scenario. 



5.3. Comparison with the early ALFALFA catalog 
release 

5.3.1. Sample selection and corrections to the line-width 

By the end of the writing of this paper, only three 
catalogs have been publicly rele ased by the ALFALFA 
collaboration, tw o in the VdR (jGiovanelli et al.1 120071 : 
iKent et all 120081) and the other one in the aVdR 
(|Saintonge et alj I2008T) . These catalogs comprise only 
~ 6% of the final volume. The first two cover an area 
from ll h 30 m < R.A. < U h and +08° < DEC < +16° 
and the last one from 22 h < R.A. < 03 h and +26° < 
DEC < +28°. Despite the limited volume, we make an 
attempt to compare our results to the observational data 
released so far. 

For such comparison, we take a sample of the AL- 
FALFA sources according to the following criteria: i) 
distances lower than 20 ft. _1 Mpc; ii) exclusion of High Ve- 
locity Clouds (HVCs) and sources with no measurement 
of HI mass; iii) removal of sources with no inclination 
measurement or with i < 30° 11 ; and finally, iv) removal 
of sources showing clear signs of interaction within a pro- 
jected angular distance smaller than the beam size of the 
Arecibo antenna (~ 3.5'). 

Criterion i) is the most stringent of all, only 14% of 
the sources in the three original catalogs fulfill it. Of the 
remaining galaxies, 81% satisfy criteria ii)-iv). The final 
sample consists of 186 galaxies in the VdR and 15 in the 
aVdR. 

It is necessary to obtain measurements for the inclina- 
tion of the sources since the 21cm line-width W50 (mea- 
sured at the 50% peak level) can be associated with V max 
only after appropriate corrections, one of them, deprojec- 
tion to an edge-on view. We discuss further below how 
we corrected the values of W§$. 

To get the inclinations, we extracted the minor-to- 
major axis ratios (b/a) by cross checking the sources with 
the GOLD Mine Database 12 (|Gavazzi et al.ll2003h . the 

Velocity measurements are not reliable for i < 30° because of 
the large deprojection correction and the uncertainty associated to 
the estimate of i. We have tested the influence of different values 
for this lower limit, 0°, 45°, 60°, and found no significant impact 
in our results. 

12 http://goldmine.mib.infn.it/ 
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Corn ell HI Archive of pointed sources 13 ([Springob et al.l 
120051 ). and the NASA/IPAC Extragalactic Database 
(NED) 14 . Using the axis ratios, we compute the incli- 
nations with the formula: 



cos 2 (i) = 



(b/a) 



i-«3 



(2) 



where go is the intrinsic axial ratio of a galaxy seen edge- 
on, we adopt qo — 0.2 as a fiducial value for all galaxies 
in our sample, in dependent of morphological type (e.g. 
iTullv et al.ll2009t ). For criterion iv) we checked, whenever 
it was possible, the optical counterpart of the sources 
using the Sloan Digital Sky Survey (SDSS) website 15 . 

The value of W50 given in the raw catalogs has al- 
ready been corrected for instrumental broadening. We 
further corrected these values for turbulent motions and 
for inclination effec t s fo llowing the procedure given in 
IVerheiien fc Sancisil (|2001l ): 



sin(i) 



■ wl 50 



1 - 2e 



2W 5 n R W, 



50,Ji"t,5Q 



' W 50,R ' 



' w fM,R 



1 



(3) 



where Wm,R is the raw value corrected for instrumental 
broadening, Wt,m = 5kms _1 and W c ,50 = 100 km s -1 . 
After these corrections, the maximum rotational veloc- 
ity of the galaxies can be estimated with reasonable accu- 
racy as W50/2. We make however the following comment 
on the estimation of V max using W50: although the HI 
line-width provides no information on the radial rota- 
tion profile of the galaxy, most of the HI gas is located 
in the outer part of disk galaxies (typically beyond three 
scale lengths) 16 . For large disks, W50 provides a measure 
of rotational velocity in the regions wh ere the rotation 
curve is alr eady flat or rising slowly (e.g. ICatinella et al.l 
2006. 120071 ). In the case of galaxies with lower velocities 
(< 75 km s _1 ), their rotation curves are typically thought 
to be still rising to the last measured point. However, 
ISwaters et al.l ()2009f ) have recently analyzed in detail a 
sample of dwarf galaxies and found that the shapes of 
their rotation curves are similar to those of more mas- 
sive galaxies; in particular, the rotation curves typically 
start to flatten at two disk scale lengths. If these re- 
sults hold for the galaxies in our sample, then we can be 
confident to use W$q to get V max without an important 
systematic underestimation. 

5.3.2. The HI velocity function 

The angular position of the final sample of galaxies 
is shown in Fig. [T^l using equatorial coordinates (upper 
and middle panels for the VdR and aVdR respectively). 
As a reference, the position of M87 is marked in the 
figure with a red star. The lower panel shows the number 
density of sources as a function of their distance to the 
MW. The prominent overdensity around 11 /i _1 Mpc is 
caused by galaxies in the vicinity of the Virgo cluster. 

13 http:/ /arecibo. tc.cornell.edu/hiarchive/ 

14 http:/ /nedwww. ipac.caltech.edu/ 

15 http://www.sdss.org 

16 We note that strongly HI deficient low mass galaxies, like 
those inside the Virgo cluster, are typically not detected by the 
ALFALFA survey. 



A similar feature is present in our CSs (see Fig. [8]). We 
note that distances in the ALFALFA catalogs take into 
ac count the peculiar v elocity field according to the model 
of lTonrv et all IpOOOh . 

The HI velocity function of the sample of galaxies in 
the VdR is shown in Fig. [T3] with red square symbols. 
The values were compute d using the T,(l/ V max ) weight- 
ing method proposed by ISchmidtl |l968). Y max is the 
volume given by the maximum distance, D max , where 
a given source could be placed and still be detected by 
the survey. This distance is determined by the sensitiv- 
ity limit of the survey, given in the ALFALFA survey 
design for a signal to n oise threshold of 6 (see eq. 5 of 
iGiovanelli etldll2005bft : 



Mpc 



0.49 



1/2 



Mhi 
10 6 M r . 



no 



200 kms" ] 



(4) 

where Mm is the HI mass associated to the source, 
7 = -1/2 for W 50 < 200 km s" 1 and 7 = -1 for 
W50 > 200kms _1 . The parameter fp quantifies the frac- 
tion of the source flux detected by the telescope's beam, 
for simplicity, we treat all sources as point sources and 
take fp = 1. Since the beam size of the Arecibo antenna 
is 3.5', the majority of the sources can be treated in this 
way. We take a fiducial value of 48 s for the integration 
time t s . For comparison with our simulations we are re- 
stricting the analysis to distances lower than 20 /i _1 Mpc, 
thus, effectively: W max = min(V( D max ), V(20/i" 1 Mpc)). 
For the majority of the sources, D max > 20h _1 Mpc, 
therefore, the volume-weights have only a minor impact 
in the number count. However, three sources in the lower 
velocity bins have a very low value of V max and their 
weights deviate strongly from the average in their re- 
spective bins; we removed them since they are not sta- 
tistically representative and can lead to a strong overes- 
timation of the velocity function. 

The results from our CSs for the corresponding field 
of view are also shown in Fig. Q2] The red dashed and 
dotted areas encompass the la regions, using Poisson 
statistics, for the predictions in the ACDM and AWDM 
cases, respectively. The VFs were constructed taking 
into account the sensitivity limit of the survey, using 
the same S(l/V maa: ) weighting method as for the ob- 
servations and renormalizing the result according to the 
fraction of galaxies that were excluded from the final ob- 
servational sample due to the 30° inclination cutoff. To 
properly apply the sensitivity limit we would need to give 
an estimate of the gas fraction, f gas — Mhi /Mdisk, in 
the modeled disks, that depends on the efficiency of gas 
transformation into stars and the gas infall history. For 
simplicity we take f gas = 1, but we note that the value 
of fgas is irrelevant for our particular analysis, see dis- 
cussion by the end of the section. The sensitivity limit 
has a relevant effect only for the low velocity disks. 

For velocities larger than ~ 80kms _1 , the VF of both 
CSs match reasonably well the observational data . This 
result is not trivial, for example, at V max — 100 km s , 
the value of the VF for halos in the whole simulated 
box is approximately an order of magnitude lower than 
the value associated to the observed sample in the VdR. 
Therefore, we confirm that the CSs are able to simu- 
late properly the overdense VdR. For velocities in the 
range 35-80 kms -1 , the ACDM simulation overpredicts 
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the value of the velocity function, increasingly for lower 
velocities. On the other hand, the AWDM simulation is 
in good agreement with the observed data, with values 
slightly higher. The low velocity end, V ma x < 35 kms -1 , 
is not suitable for comparison since both, simulations 
and observations are not complete for the lower veloci- 
ties. The minimum halo mass in our simulations that is 
reliable, according to a comparison of the MF with the- 
oretical expectations in the low mass end (see Fig. 2]), 
sets a lower limit of completeness for V rnax . For the 
ACDM simulation, this mass is ~ 10 9 /i _1 M0, corre- 
sponding to Vmax ~ 24 kms -1 ; in the AWDM case is 
given by the limiting mass 3 x 10 9 ft. -1 M Q related to 
Vmax ~ 29 kms -1 . The typical sensitivity limit for the 
ALFALFA survey goes down to M HI = 10 7 ft -1 M Q for 
distances up to 20 /i -1 Mpc, see eq. @, a source with 
this mass can have a V max value within a broad range, 
due to the natural scatter on the V max — Mm relation. 
An examination of our galaxy sample shows that this 
range is ~ 13 — 35 kms -1 , which puts a lower limit of 
~ 35 kms -1 for the completeness of the sample. Tak- 
ing the theoretical and observational limits into account, 
a consistent comparison can only be done for values of 
Vmax above ~ 35 kms -1 . 

So far we have used a constant disk baryon fraction 
fdisk — 0.03 to populate the dark matter halos with 
galaxy disks. A more realistic approach would be to 
include the effects of SN feedback in fdisk- This in- 
clusion would only have an important effect for low ve- 
locity halos if the gas outflow produced by SN is suf- 
ficient to deplete the forming galaxy of gas and push 
it below the sensitivity limit of ALFALFA. We now in- 
corporate thjs_effie£^jas2n£_ttie_disk galaxy evolutionary 
model of iDutton fc van den Bosch] (|2008h for the case- 
where the outflow of gas is energy driven: the kinetic 
energy of the wind is a fraction €efb = 0.25 of the 
kinetic energy produced by SN. In this case, fdisk be- 
comes a str ong function of halo mass (see upper panel 
of Fig. 8 in IDutton fc van den Bosch! 120081 ). We adopt 
the fitting function obtained by the authors for the en- 
ergy driven feedback model (see their eq. 43 and table 3) 
and apply it to our model. The results appear in Fig. [13] 
as dashed and dotted lines for the ACDM and AWDM 
CSs. The net impact is a marginal reduction of the VF 
at the low velocity end. A more drastic reduction would 
be possible with a model that significantly lowers the 
value of fgas, which we have set to one for simplicity. 
For instance, the lowest value we have set for compari- 
son, Vmax = 35 kms -1 , is related on average to a halo 
of 6 x 10 9 ft -1 M , which in the case of the SN feedback 
model corresponds to fdisk = 0.02, thus the modeled disk 
would lie below the sensitivity limit of the survey only if 
fgas is lower than 0.08. In our scheme the gas fraction 
only determines if a given halo should be counted or not, 
thus, the previous result indicates that the specific value 
of fgas bas no impact unless the galaxies we are aiming 
to compare with have very low gas fractions. HI defi- 
cient low mass galaxies are typically present inside galaxy 
clusters like Virgo and are unlikely to be detected by AL- 
FALFA; they also have no counterpart in our simulations 
since we have dealt exclusively with halos, excluding the 
subhalos within them. Galaxies in the field are less likely 
to have such low gas fractions. Nevertheless, we explic- 
itly checked that this is indeed the case for the sample 



of galaxies in the VdR. For that purpose, we identified 
the SDSS optical counte rparts of the HI sources an d used 
the method described in lBlanton fc Roweii (|2007l ) to es- 
timate the stellar masses of 73% of the galaxies in the 
final VdR sample. We found that the majority of these 
galaxies have large gas fractions, approximately 90% of 
the ones with V ma x < 100 kms -1 have f gas > 0.1. 

In Fig. [14] we show analogous results for the aVdR. Al- 
though in this case the statistical significance of the ob- 
served sample is much lower (only 15 sources), the main 
results observed in the case of the VdR are reproduced, 
indicating that simulations are able to probe properly 
this region as well. 

6. DISCUSSION AND CONCLUSIONS 

N-body simulations with constrained initial conditions, 
set-up to resemble the spatial distribution of dark mat- 
ter in the local Universe, are a powerful tool to make de- 
tailed comparisons between predictions of the dark mat- 
ter paradigm and observational evidence in our local en- 
vironment. 

Using the algorithm of constra ined realizations devel- 
oped by [H offman & Ribak (1991 ) and the approach laid 
out in lKravtsov et alJ (|2002f ) and lKlvpin et al.) <|2003f l. we 
have run a pair of CSs that incorporate nearby observa- 
tional data sets as input for their initial conditions. One 
simulation follows the evolution of structure in a CDM 
model and the other one is based on a thermal WDM 
particle with a mass of IkeV resulting in an effective fil- 
tering mass of ~ 10 10 h~ 1 M Q . 

After an appropriate choice of the coordinate system, 
the simulations are able to reproduce the overall spa- 
tial distribution of the most significant structures within 
20 /i -1 Mpc of the Local Group, namely, the LSC, includ- 
ing the Virgo cluster, as well as the Fornax cluster laying 
in the opposite direction (see section 3). 

The mass and velocity functions of halos in the whole 
simulated boxes are consistent with theoretical expec- 
tations. The ACDM case follows closely the estimates 
from the Sheth & Tormen formalism, except in the high- 
est mass end due to the influence of the LSC. For the 
AWDM cosmogony, the results lie close to the ACDM 
case for halos with masses higher than the filtering mass; 
for lower masses, the mass and velocity functions flatten 
and then rise due to spurious numerical fragmentation 
for masses < 3 x 1 9 fe~ 1 M f 7) confirming t he limiting mass 
formula given bv I Wang fc Whltel (p007h . The mass res- 
olution of the CSs, rriDM = 1-63 x 10 7 h~ 1 M Q , allows us 
to derive robust results for halos with masses larger than 
this limiting mass, corresponding to maximum rotation 
velocities of 24 kms -1 . 

For the local environment, the general prediction of 
the CSs is that the region within 20 /i -1 Mpc of the LG 
is overdense by a factor of ~ 2 compared to the MF 
of the whole simulated boxes (see Fig. [6]). Such result 
goes in the same direction as a recent result reported by 
iTikhonov fc Klypinl (|2009f ). showing that the luminosity 
function in a sample of galaxies within 8 h~ 1 Mpc is larger 
than the universal luminosity function by a factor of 1.4. 

Since we have shown that our CSs are capable of repro- 
ducing the abundance of halos in the local environment, 
we have obtained predictions for the VF of halos in the 
field of view, within 20 /i -1 Mpc, that is being surveyed 
by ALFALFA. The VF has the important advantage over 
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the MF that it can be compared more directly with ob- 
servational data since it avoids the problem of relating 
halo masses to galaxy luminosities. 

After completion, the ALFALFA survey will detect HI 
sources covering 7000deg 2 of the sky in two different re- 
gions. In this work, we have referred to these regions, 
additionally constrained to 20/i _1 Mpc, as the Virgo- 
direction region, VdR, and the anti- Virgo-direction re- 
gion, aVdR. Our CSs predict that the VF of halos in the 
VdR exceeds the universal velocity function by a factor 
of ~ 3. The VF in the aVdR is only slightly underdense 
(~ 10%) than the universal value (see Fig. [9]). 

We have used a simplified model to populate our halos 
with disk galaxies. It only incorporates the dynamical 
effects of the disks. With this model, we are able to pre- 
dict the VF of disk galaxies in the two regions explored 
by ALFALFA. Although the survey is not complete yet 
we have compared our predictions and the results from 
a sample of galaxies taken from the catalogs released so 
far, which cover only 6% of the total planned volume. 

For velocities in the range between 80 kms -1 and 
300km s _x , the VFs predicted for the ACDM and A WDM 
simulations agree quite well with the VF of the sample of 
galaxies. For the VdR, this result is particularly encour- 
aging and reassures the confidence in our CSs to properly 
simulate the local environment: despite the small volume 
used for comparison, the simulations are able to predict 
the shape and normalization of the VF in the high ve- 
locity regime. In the VdR, the normalization is an order 
of magnitude larger than for the universal VF. 

For velocities larger than the minimum mass we can 
trust for comparison, ~ 35 kms -1 , and lower than 
80 kms -1 , the predictions agree well for the AWDM 
cosmogony (contrar y to recent claims, see section 9 of 
iBlanton et al.l [20081 ). with the VF being approximately 
flat in this regime. The ACDM model, however, predicts 
a steep rise in the velocity function towards low veloci- 
ties; for V m ax ~ 35 kms -1 , it predicts ~ 10 times more 
sources than the ones observed. Using the same set of 
simulations Tikhonov et al. (2009), in preparation, found 
that the observed spectrum of mini- voids in the local vol- 
ume is in good agreement with the AWDM model but 
can hardly be explained within the ACDM scenario. 

Although we have only explored a simplified model to 
populate our halos with disk galaxies, our results indicate 
a potential problem of the ACDM paradigm in the low- 
velocity regime of dwarf galaxies. Nevertheless, there are 
several issues that need to be addressed before reaching 
a strong conclusion. 

On the observational side, the sample of galaxies we 
have analyzed comprises only 186 galaxies in the VdR. 
It is of key importance to obtain the VF for the complete 
volume of ALFALFA and check if the flattening at low ve- 
locities is reproduced. Another important issue is to de- 
termine to which extent the value of W50 is representative 
of the maximum rotational ve locity for the least m assive 
galaxies. A recent analysis bv lSwaters et alJ (2009) indi- 
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cates that the rotation curves for late-type dwarf galaxies 
are similar to those of more massive galaxies, starting to 
flatten at two disk scale lengths, thus, W50 would not un- 
derestimate strongly the value of V ma x for these galaxies 
since the HI gas typically extends beyond three disk scale 
lengths. 

From a theoretical perspective, astrophysical phenom- 
ena such as SN feedback and UV photoionization play 
an important role to deplete gas from low mass ha- 
los. The relevant influence of these effects on the VF 
comes from the typical sensitivity limit of the survey, 
allowing to detect HI masses larger than 10 7 h~ 1 M@ 
at distances D < 20ft~ 1 Mpc. We have explored the 
fo rmer of these effects by using the model presented 
in iDutton fc van den Bo sch l|2008l) and found little im- 
pact on th e VF. Regarding UV background heating, 
iHoeft et al.l (|2006f ) found that this mechanism is not very 
efficient in evaporating all baryons from dwarf sized ha- 
los below a characteristic mass scale of 6 x 10 9 ft _1 MQ 
(V max ~ 35 kms -1 ). Since according to our findings, a 
strong suppression of gas is needed for halos with charac- 
teristic velocities V max < 60km s -1 to flatten the velocity 
function in the ACDM case, it is unlikely that UV pho- 
toionization could account for it. 
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Fig. 1. — Unconstrained linear power spectrum at z=0 for the ACDM and AWDM (mwDM = 1 kcV) cosmologies (blue and red 
respectively). The associated Nyquist frequency and fundamental mode for our set of simulations are represented with vertical lines to the 
right and left in the figure, respectively. 



TABLE 1 

Main properties of the LG, Virgo and Fornax. The first line 

IS THE MASS OF THE OBJECT AND THE SECOND THE RELATIVE 
DISTANCE TO THE LG. DISTANCES ARE GIVEN IN /l _1 MPC AND 
MASSES IN lO 12 h~ 1 M0. 



Object/Case 


LG(MW+M31) a 


Virgo 


Fornax 


ACDM 


1 


100 


27 




0.9 


10.8 


10.5 


AWDM 


1 


94 


31 




0.5 


10.2 


10.6 


Obs. 6 


1.6 C 


204 


31 




0.6 


11.4 


15.0 



"The distance is in this case the relative distance between MW and 
M31 

'Cluster values according to Girardi ct al. (1998) 
^Cen tral values of the ma sses estimates acco r ding to IXue et aLl 
| |2008T) for the MW and from IEvans~ Wilkinson (200Cj) for M31 
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Fig. 2. — Left panel: projected dark matter distribution for the ACDM constrained simulation at z = 0. The slice has a thickness 
of 8/i _1 Mpc and it is centered at Z = 24h _1 Mpc (in Box coordinates). It encompasses part of the supergalactic plane. The Local 
Superclustcr is the filament crossing the image horizontally. The locations of the Local Group and the Virgo cluster are marked in the 
figure. Right panel: same as on the left for the AWDM simulation. 




MW-like halo in the ACDM simulation for coordinate system SGzero- The color scale represents a smoothed mass-weighted count of the 
number of halos per pixel, going from blue to red for lower to higher counts. The angular positions of all halos with masses larger than 
5 X 10 9 /i — 1 Mq are displayed with black dots. The boxes marked with solid black lines in the center and in the sides are the VdR and 
aVdR respectively. The locations of the simulated Virgo and Fornax clusters appear as black circles, the correspondent locations of the 
real clusters appear as black squares. Right panel: same as on the left for coordinate system SG m i„. 
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Fig. 4. — Upper panel: differential MFs for the ACDM and AWDM constrained simulations (blue and red solid lines respectively). 
Analytical predictions following the Sheth & Tormen formalism for both cosmologies are shown as dashed lines. Lower panel: ratio of the 
MFs to the value of the MF given by the Sheth & Tormen formalism for the ACDM cosmology. The solid and dotted vertical lines mark 
the values of the filtering and limiting mass of the AWDM simulation: 1.1 X 10 10 h~ 'Mq and 3 X 10 9 /i -1 Mq respectively. 
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Fig. 5. — The same as Fig. [4] but for the VFs instead of the MFs. In this case no analytical prediction is shown for the A WDM cosmology. 
The vertical lines indicate the maximum velocities corresponding to the filtering and limiting masses, with values of 36kms _1 and 24km s — 1 
respectively. 
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Fig. 6. — Ratio of the differential MFs for the ACDM (upper panel) and AWDM (lower panel) CSs to the prediction of the S-T formalism 
for the ACDM cosmology. The solid lines represent the result for the whole simulation boxes (see Fig. [3J. The blue dashed lines are for 
the S-T prediction. The other lines are for spheres of different radii centered at the LG. The curve for the 20/i -1 Mpc radius sphere is 
highlighted with a thick dash-dotted line. 
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Fig. 7. — Projections of the halo distributions in the ALFALFA field of view onto the planes of the supergalactic coordinate system 
SG m i n . Only halos within 20 h~ 1 Mpc from the MW halo and with masses larger than 5 X 10 9 /i — 1 Mq are shown. The red and blue dots 
represent halos inside the VdR and aVdR respectively. The left and right panels are for the ACDM and AWDM simulations, respectively. 
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Fig. 8. — Upper panel: number density of halos, in the ACDM simulation, as a function of distance to the LG for the halos displayed 
in Fig. [7] The values are normalized to the number density within the 20 h~ 1 Mpc sphere. Red and blue lines are for the VdR and aVdR 
respectively. The black line gives the radial number density profile for the whole sphere. Lower panel: same as on the left for the AWDM 
simulation. 
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Fig. 9. — Upper panel: differential VFs for the ACDM (solid lines) and AWDM (dotted lines) constrained simulations in the VdR (red) 
and aVdR (blue). The black (solid and dotted) curves are for the sphere of 20/i -1 Mpc radius centered at the LG. The S-T prediction, 
only for ACDM cosmology, is shown as a black dashed curve. Lower panel: ratio of the VFs for the simulations to the VF given by the 
S-T formalism for the ACDM cosmology (same line styles as in the upper panel). The solid and dotted vertical lines mark the values of 
the filtering and limiting velocities of the AWDM simulation as in previous figures. 
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Fig. 10. — Ratio of the differential VF of the ACDM simulation for the VdR to the VF in a sphere of 20/i _1 Mpc radius. The red 
zone encompasses the value of these ratios for different rotations of the supergalactic coordinate system. The results for two of these 
rotations, corresponding to the coordinates systems SGzero and SG m i n denned in section[3] are highlighted as black dashed and solid lines 
respectively. Also shown are the ratios in the aVdR for these two rotations (blue lines). 
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Fig. 12. — Upper panel: distribution in equatorial coordinates of the sample of HI sources taken from the ALFALFA public catalogs 
released so far in the VdR. The position of M87 in the Virgo cluster is marked as a red star. Middle panel: sample of sources in the aVdR. 
Lower panel: the number density of sources as a function of their distance to the MW (red and blue for the VdR and aVdR respectively). 
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Fig. 13. — VF for the sample of galaxies in the VdR taken taken from the ALFALFA catalogs (square symbols with error bars). 
Predictions from our CSs, using fdisk = 0.03, for the observed field of view appear as the dashed (ACDM) and dotted (AWDM) red 
areas, delimited by Poisson error bars. The dashed and dotted lines are predictions using a model where fdisk is a function of halo mass 
incorporating the effects of SN feedback. The sensitivity limit of the survey has been included in the results. The vertical solid line marks 
the value of V m ax down to which the simulations and observations are both complete. 
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Fig. 14. — The same as Fig. 1131 but for the limited survey coverage of the aVdR. 



